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SENSITIVITY  ANALYSIS  OF  THE  SAGNAC-EFFECT 
OPTICAL-FIBER  RING  INTERFEROMETER 

INTRODUCTION 

Recently,  R.B.  Brown  [1],  in  a study  of  inertial  rate  sensing,  suggested  the  use  of  a 
Sagnac-effect,  fiber,  ring  interferometer.  Through  the  use  of  a multiple-turn  fiber  loop,  high 
sensitivities  appeared  to  be  feasible  if  fiber  optic  attenuation  could  be  reduced.  Since  that  time 
(1968),  the  development  of  single-mode,  extremely  low-loss  optical  fibers  [2]  has  progressed 
rapidly,  making  fiber-optic,  Sagnac  interferometers  possible.  Vali  and  Shorthill  [3,4]  were  the 
first  to  demonstrate  single-fiber  interferometer  gyro  for  rotational  rate  sensing.  Fiber  inter- 
ferometers using  two  fibers  represent  a second  class  of  devices  that  appear  to  be  suitable  as 
highly  sensitive  pressure  [5],  stress,  and  temperature  sensors.  In  this  report,  an  analysis 
describing  single  fiber  interferometer  gyros  is  presented  that  provides  the  theoretical  senstivities 
of  these  devices  in  practical  situations. 

In  their  first  paper,  Vali  and  Shorthill  [3]  gave  a preliminary  estimate  of  fiber  interferome- 
ter sensitivity;  however,  their  estimate  was  not  intended  to  completely  describe  the  perfor- 
mance of  these  devices  and  consequently  did  not  accurately  describe  the  experimental 
configurations  most  commonly  encountered.  Specifically,  they  did  not  take  into  account  the  ubi- 
quitous scattered  light  in  the  long  optical  fiber  as  one  of  the  noise  sources.  Furthermore,  the 
formula  they  used  to  estimate  the  photon-noise-limited  sensitivity  of  the  interferometer  was  not 
derived  for  the  experimental  optical  configurations  depicted  in  their  papers;  [3,4],  but  rather, 
was  directly  taken  from  a paper  by  G.E.  Moss,  et  al.  [6].  Moss,  et  al.  considered  a heterodyne 
ac  phase-detection  technique,  whereas  the  optical  detection  used  by  Vali  and  Shorthill  [3,4], 
was  a dc  (or  homodyne)  technique.*  The  results  for  an  ac  detection  scheme  do  not  necessarily 
apply  to  dc  schemes  because  one  has  to  consider,  in  each  case,  different  noise  sources  that 
affect  the  outcome  of  a sensitivity  analysis,  and  hence  its  conclusions. 

This  report  presents  some  results  of  a sensitivity  analysis  of  Sagnac-effect,  optical-fiber, 
ring  interferometers  with  detailed  considerations  of  signals  and  noises.  For  simplicity  as  well  as 
to  provide  a focus  of  discussion,  the  experimental  optical  fiber  configuration  proposed  by  Vali 
and  Shorthill  was  analyzed  and  compared.  In  the  following  sections,  we  will  first  derive  per- 
tinent formulas  for  signal-to-noise  ratios,  and  then  discuss  noise  components  due  to  Rayleigh, 
Brillouin,  Mie,  or  core-cladding  interface  light  scattering.  These  will  be  followed  by  a discus- 
sion of  particular  problems  and  promising  improvements.  Finally,  the  limiting  sensitivities  at 
both  low  and  high  input  powers  are  given. 


Manuscript  submitted  June  8,  I97S.  Resubmitted  July  3,  1978. 

•The  modulation  frequency  of  this  interferometer  was  set  as  high  as  20  kHz  so  that  low-frequency  noise  effects  such  as 
mechanical  vibration-induced  noise  (l//noise  in  detection)  were  eliminated  by  narrow  bandpassins.  * L.  Forward  has 
pointed  out  the  inadequacy  of  using  the  result  of  Ref.  6 in  a dc  detection  system  (private  communication). 
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The  development  of  a ring  interferometer  gyroscope  that  is  lightweight,  low-cost,  and 
compact  offers  an  alternative  to  the  nearly  maturing  ring  laser  gyros  and  well-matured  mechani- 
cal gyros,  provided  that  stable,  phase  biasing  techniques  can  be  developed.  These  biasing  tech- 
niques described  in  this  report  are  required  so  that  the  gyro  can  be  operated  in  the  region  of 
maximum  sensitivity. 


SIGNAL-TO-NOISE  ANALYSIS 
Operational  Description 

The  optical  configuration  on  which  the  present  analysis  is  based  is  shown  in  Fig.  1.  The 
optical  source  is  a single-mode,  stabilized,  coherent  (gas  or  semiconductor)  laser  with  output 
power  Pq,  and  wavelengths  in  the  visible  or  near-infrared  region.  The  laser  output  beam  is 
assumed  to  be  well  collimated,  with  uniform  phase.  As  the  beam  first  passes  through  a lossless 
beam  splitter  BS|,  which  has  a power  reflection  coefficient  a{,  its  power  is  attenuated.  The 
transmitted  part,  /Vl-ai),  proceeds  through  a lossless,  nondispersive  medium  and  reaches  the 
second  beam  splitter  BS2,  which  is  designed  to  give  a SO/SO  split  of  the  power  at  an  incident 
angle  of  45°:  (1/2)  Pod  — a ,)  is  the  power  of  the  reflected  as  well  as  that  of  the  transmitted 
beam.  Each  beam  is  then  launched  into  one  end  of  the  single-mode  optical-fiber  coil  of  radius 
R by  a focusing  lens  system,  such  as  a microscope  objective,  with  the  optimal  //number  to 
ensure  maximum  power-coupling  efficiency  C.  In  continuous-wave  operation,  there  will  be  two 
waves  propagating  in  the  whole  length  L of  the  fiber  simultaneously  but  in  the  opposite  direc- 
tions. If  the  total  loss  in  power  through  the  optical  fiber  follows  the  exponential  law  with  an 
attenuation  coefficient  ar(dB/km),  and  if  the  loss  is  reciprocal,  i.e.,  if  the  loss  is  the  same 
irrespective  of  the  wave  propagation  direction,  then  the  power  of  each  emerging  beam  will  be 
(l/2)/*0C(l  -a,)  e "r' , which  we  will  assume  to  be  100%  intercepted  and  recollimated  by  the 
other  focusing  lens  system.  The  beam  in  the  clockwise  direction  (CW)  is  to  undergo  one  more 
reflection  and  the  beam  in  the  counterclockwise  direction  (CCW),  one  more  transmission, 
before  they  combine  to  form  a fringe  pattern  F t.  The  power  for  each  beam  is  then 
(l/4)C/»0(l-a,)e"“ri. 
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As  they  continue  toward  the  BS,,  the  remaining  halves  of  the  beams  are  reflected  from  it 
to  form  the  second  fringe  pattern  £2.  Their  powers  will  be  (l/4)C/,ot»|(l -a,)  "T!  This 
fringe  pattern,  formed  by  waves  that  have  undergone  different  numbers  and  orders  of 
reflections  and  transmissions  from  those  of  F\ , has  a phase,  in  general,  different  from  that  of 
F\.  If  the  difference  is  180°,  the  two  fringe  patterns  are  said  to  be  complementary.  In  what 
follows,  we  will  assume  that  this  is  the  case. 

A well-adjusted  optical  system  such  as  this,  if  in  inertial  motion,  will  give  rise  to  two 
infinite-width  fringe  patterns  that  have  uniform  but  different  brightnesses.  However,  as  soon  as 
the  system  is  set  into  rotation,  concentric  interference  rings  are  formed  due  to  the  phase 
difference  between  the  CW  and  the  CCW  waves  because  of  the  Sagnac  effect  [71.  When  the 
rotational  axis  is  parallel  to  the  axis  of  the  optical  fiber  coil  of  N turns  and  of  constant  loop  area 
A,  the  phase  difference 
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where  X is  the  free-space*  optical  wavelength,  c is  the  free-space  speed  of  light,  and  11  is  the 
rate  of  rotation.  For  a circular  coil,  A *■  n R2  and  L — 7ttR.N\  therefore, 
A4  - [4  n LR/Occ) ) Cl.1  The  quantity  in  the  square  brackets  is  the  sensitivity  of  the  phase 
difference  to  the  rotational  rate.  It  is  clear  in  principle  that,  to  increase  the  sensitivity,  we 
should  increase  the  total  length  and  radius  of  the  coil,  and  use  as  short  an  optical  wavelength  as 
possible.  However,  in  practice,  packaging  criteria  limit  the  size  of  R\  optical  fiber  loss  sets  an 
upper  bound  on  the  length  L;  nonlinear  damage  effects  forbid  the  use  of  very  high  power  laser 
sources;  and  the  signal  strength,  the  scattered  light,  and  the  quantum  efficiency  of  the  detector 
ultimately  limit  the  system  sensitivity.  In  other  words,  system  tradeoffs  must  be  studied.  The 
sensitivity  analysis  done  below  is  thus  central  to  the  trade-off  study. 

The  Signal 

The  ring  interferometer  described  above  is  not  sufficient  for  rotational  rate  sensing 
because  one  has  to  extract  the  phase  information  from  the  fringe  patterns.  Although  fringe 
counting  techniques  [8,9]  have  been  used  to  extract  the  phase  by  spatial  sampling  of  the  fringe 
pattern,  they  involve  either  a modification  of  the  interferometer  setup  or  a multisensor  fringe 
locator  with  rather  elaborate  electronic  logic  circuitry.  A simpler  technique,  as  shown  in  Fig.  2, 
is  to  use  a single  photodetector  focused  on  the  center  of  the  fringe  pattern  such  that  the  area  in 
the  fringe  pattern  in  the  field-of-view  of  the  detector  is  small  compared  with  the  width  of  a 
fringe  to  ensure  uniform  optical  intensity  across  the  detector  surface. 

The  fringe  intensity  depends  on  the  interference  (i.e.,  coherence)  among  the  optical  fields 
emerging  from  one  end  of  the  fiber  and  those  from  the  other  end.  Let  Elcw  be  the  sum  of  all 
possible  optical  fields  from  one  end  in  the  CW  direction;  then,  in  general,  £,cw  - Ecoh  + 
£jncoh,  where  £coh  represents  the  sum  of  coherent  optical  fields,  and  £mcoh  represents  that  of 


•By  "free  space"  here,  we  mean  either  a vacuum  or  a homogeneous  nondispersive  medium  in  which  the  optical  system 
is  submerged 

t The  fringe  shift  in  a fringe  pattern  caused  by  the  phase  difference  is  just  \<t>/2n. 
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PHOTODETECTOR 

Fig.  2 — Optical  mixing  technique  for  the  phase 
difference  measurement 


the  incoherent  fields.  Similarly,  £tccw  — £coh  + £incoh,  in  the  CCW  direction.  Then  the  instan- 
taneous intensity  in  the  fringe  is  given  by 

(^coh  + ^incoh)  (^coh  4"  ^incoh) 


I 


C’ 

lew  X-tCCW 

“ ^coh^coh  £ 


incohccoh 


^"coh^incoh  "h  ^inenh^inc 


incohcincoh 


(1) 


If  this  instantaneous  intensity  is  allowed  to  fall  onto  a detector,  it  will  generate  photoelectrons, 
which  in  the  end  produce  a current, 

/ - Dl,  (2) 

where  D is  the  detector  conversion  factor.  In  Eq.  (1)  or  (2),  after  we  take  the  time  average 
over  an  interval  large  compared  with  the  period  T — 2ir/cu,  we  will  get  terms  of  autocorrela- 
tions (i.e.  intensities)  of  both  coherent  and  incoherent  components,  as  well  as  terms  involving 
crosscorrelations  between  coherent  components.  The  crosscorrelations  among  the  incoherent 
components,  or  between  coherent  and  incoherent  components,  will  vanish  by  definition.  Only 
the  crosscorrelations  between  coherent  components  carry  definite  phase  information.  Therefore, 
the  photocurrent  /,  after  time  averaging,  will  consist  of  a dc  term  involving  all  the  intensity 
terms  and  an  ac  term  involving  sinusoidal  functions  of  phases  between  interfering  coherent 
beams.  (The  specialization  of  this  to  a detector  with  finite  response  time  in  a ring  interferome- 
ter will  be  discussed  below.)  The  photocurrent  is  then  amplified,  bandpassed,  and  integrated. 
This  is  the  optical  mixing  technique  [10]  and  is  well  studied  in  areas  of  laser  communication 
systems  [11].  We  note  that  this  optical  configuration  does  not  contain  any  external  modulation. 
Therefore,  the  detection  is  direct  current  in  nature.  It  differs  from  the  conventional  homodyne 
detection  scheme,  in  that  there  is  no  adjustable  "reference  beam."  The  beats  are  formed  by  sig- 
nal beams  themselves.  Because  of  these  special  features  in  the  proposed  optical  configuration, 
we  cannot  read  out  the  phase  difference  directly. 

The  optical  mixing  output  signal  is  formed  by  the  two  directly  transmitted  beams.  The 
amplitude  of  the  signal  depends  strongly  on  the  coherence  (spatial  and  temporal)  and  polariza- 
tion states  of  the  two  beams.  These  states,  in  turn,  are  influenced  by  scattered  light,  optical 
misalignment  due  to  mechanical  vibrations,  nonuniform  temperature  and  stress  fields  in  the 
optical  fiber,  etc.  To  simplify  the  analysis,  we  assume  that  the  optical  components  are  rigidly 
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tied  down,  the  temperature  and  stress  fields  are  homogeneous  and  stationary,  the  spatial  and 
temporal  coherence  of  the  beams  are  not  degraded  when  transmitting  through  'he  single-mode 
optical  fiber,*  and  that  polarization  vectors  of  the  two  beams  remain  aligned.!  We  may,  as 
alternative,  consider  that  the  two  beams  are  unpolarized.  We  also  assume  that  the  scattered- 
light  components  are  completely  uncorrelated  with  the  directly  transmitted  beams  so  that  the 
only  effect  to  be  considered  here  is  their  contribution  to  the  background  noise. 


Noise  Sources 

To  understand  the  noise  sources  that  are  possibly  present,  we  should  consider  the  follow- 
ing: 

1.  Scattered  light  in  the  optical  fiber 

2.  Shot  noise,  background  noise,  thermal  noise,  generation-recombination  noise,  flicker  or 
1//  noise,  dark  current  noise,  and  amplifier  noise 

3.  Fluctuations  in  the  laser  light 

4.  Environmental  noise  sources  such  as  thermal  variations,  acoustic  agitations,  and 
mechanical  vibrations  or  any  other  large-scale  perturbations  such  as  earth’s  magnetic  field  varia- 
tions, etc. 

Sources  (1)  and  (2)  are  intrinsic  because  they  are  not  related  to  any  environmental  fac- 
tors. The  fluctuations  in  laser  light  will  cause  amplitude  modulation  in  the  signal;  but,  we  will 
see  below  that  using  a differential  scheme  involving  the  two  signals  coming  from  two  detectors 
set  in  the  two  fringe  patterns  can  eliminate  the  common  mode  of  the  fluctuations.  The  scheme 
is  also  effective  in  eliminating  all  the  common-mode  environmental  noises.  This  is  the  salient 
feature  of  the  optical  configuration  under  consideration.  Hence,  in  what  follows,  we  will  con- 
centrate only  on  the  intrinsic  noise  sources. 

The  scattered  light  in  the  optical  fiber  consists  of  various  components  caused  by  different 
scattering  mechanisms.  Parts  of  the  scattered  light  are  trapped  in  the  core  and  guided  to  both 
ends  of  the  fiber  and  contribute  to  the  noises;  other  parts  are  scattered  out  of  the  optical  fiber 
and  lost  forever.  The  major  linear  mechanisms  are  Rayleigh  scattering,  Mie  scattering,  core— 
cladding  interface  scattering,  and  Brillouin  scattering. t Of  all  these  components,  only  the  Bril- 
louin  component  has  a small  frequency  shift  at  the  backscattering  direction;  the  rest  are  oscillat- 
ing at  the  same  frequency  as  the  input  laser  beam.  Whereas  the  intensities  of  the  components 
due  to  Rayleigh  and  Brillouin  scattering  are  symmetrically  distributed  with  respect  to  the  for- 
ward and  backward  directions,  those  due  to  Mie  and  core— cladding  interface  scattering  are 

•This  is  certainly  not  true  for  a multimode  fiber,  for  it  has  been  observed  that  speckled  patterns  are  formed  at  the  out- 
put of  such  a fiber 

fFor  linearly  polarized  input  waves,  the  output  waves  from  a long,  single-mode  fiber  may  acquire  some  degree  of  ellip- 
licily  in  their  polarization  states  due  to  stress  birefringence  or  deviation  of  the  core  cross  section  from  a perfectly  circu- 
lar shape  However,  these  effects  are  most  likely  reciprocal  Therefore,  in  keeping  with  the  assumption  of  reciprocity 
in  liber  losses,  we  feel  that  it  is  reasonable  at  this  lime  to  put  aside  the  question  on  polarization 

(There  is  Raman  scattering,  which  has  a much  larger  frequency  shift  Because  its  scattering  coefficient  is  10  to  20  limes 
smaller  than  Brillouin  scattering,  we  will  not  consider  it  here. 
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mainly  one-sided,  i.e.,  only  in  the  forward  direction.  Therefore,  in  the  CW  direction,  we  will 
have,  aside  from  the  direct  CW  beam  (its  field  denoted  by  £cw),  a Rayleigh  forward  component 
£RI.,  a Brillouin  forward  component  £B0,  and  a strong  forward-peaked  component  (£FP),  which 
are  associated  with  CW  beam,  and  furthermore,  a Rayleigh  backward  component  £RB  and  a 
Brillouin  backward  component  £B„,  which  are  induced  by  the  CCW  beam  traveling  in  the  oppo- 
site direction.  Thus,  we  have  a total  of  six  fields  propagating  in  the  same  one  direction.  Simi- 
larly, we  have  another  total  of  six  in  the  other  direction.  These  are  illustrated  in  Fig.  3. 


CW  ► 

v-  — i 

* rrui 
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FREQUENCY 

^ FIBER 

PHASE  FIELD  SYMBOL 

(1) 
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*CW  i*CCW 
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FORWARD  PEAK 

u 

RANDOM 

E FP 

(3) 

BRILLOUIN  SCATTERING 
FORWARD 

u 

RANDOM 

E 80 

(4) 

BACKWARD 

U) 

RANDOM 

Ebw 

(5) 

RAYLEIGH  SCATTERING 
FORWARD 

tu  *A  ft  q 

RANDOM 

Erf 

(6) 

BACKWARD 

u> 

RANDOM 

erb 

THE  OTHER  SIX  FOR  THE  CCW  CASE 


Fig.  3 — Direct  optical  beams  in  clockwise  (CW)  and  counterclockwise  (CCW) 
directions  in  a fiber  inducing  other  scattered  components:  forward  peak,  Bril- 
louin, and  Rayleigh  Excepi  for  the  direct  beams,  the  phases  of  the  others  are 
random. 


As  these  scattered-light  components  are  inevitable  in  the  optical  wave  transmission 
through  any  scatter-dominated  fibers,  and  as  they  are  always  encountered  in  complete  systems 
noise  analysis,  we  shall  devote,  in  a separate  section,  a more  detailed  discussion  to  their  origin, 
characteristics,  and  distribution  laws,  as  well  as  the  magnitudes  of  their  scattering  coefficients. 

Signal-to-Nolse  Ratio 

The  six  optical  fields  in  the  CW  direction  will  mix  with  the  other  six  in  the  CCW  direction 
on  the  surface  of  a photodetector  placed  at  the  center  of  the  fringe  pattern  (say,  £().  Since  a 
photodetector  is  a square-law  detector,  therefore,  there  are  36  terms  including  all  the  squared 
and  cross-product  terms.  But,  in  steady-state,  continuous-wave  operations,  all  the  scattered 
component  fields  will  be  incoherent  (or  of  random  phases)  because  each  is  composed  of  contri- 
butions from  randomly  distributed  scattering  centers  throughout  the  optical  fiber.  Since  the 
Brillouin  backscattered  beams  have  a frequency  shift  of  the  order  of  23  GHz,  beating  between 
them  and  the  direct  beams  may  occur.  But,  a beat  of  25  GHz  is  higher  than  any  realistic  detec- 
tor response  [12],  and  it  will  not  be  detected.  Therefore,  the  final  detector  current  output,  after 
time-averaging  over  its  time  constant  r (where  r » T - 2 ir/<o),  as  given  by 


I 
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7,  - — f /' j (/ ) dt  (3) 

J <*  o 

(where  the  overbar  indicates  time  average  and  the  subscript  1 means  photocurrent  obtained  in 
fringe  pattern  £,),  will  contain  only  one  cross-term  involving  the  two  coherent  signal  beams. 
However,  the  squared  terms  of  all  six  fields  survive  the  averaging  operation  and  show  up  in  the 
dc  term  of  the  detector  current.  If  the  laser  fields  at  the  output  ends  of  the  fiber  are 
represented  by  sinusoidal  functions  £cw  — ACVI  cos  (o>r— </>cw)  and  £ccw  * Acc„  cos  (cuf— <£ctw), 
then 

£cw  ” y A cw  A ccw  cos  (tt  ~t~  A0) 

— — — A cw  A ccw  cos  (A<£),  (4) 

where 

A cw  and  A ccw  are  the  amplitudes 
to  is  the  laser  circular  frequency 
i is  time 

<f>cw  and  </>ccw  are,  respectively,  the  phases  associated  with  the  direct  CW  and 
CCW  traveling  waves,  and  A<£  - <£cw  - 0ccW 
7T  is  a phase  shift  due  to  reflection  from  the  beam-splitter  BS2,  provided  that  it 
has  a dielectric  surface,  and  the  numerical  factor  1/2  is  due  to  the  power 
reduction  upon  reflection  from  or  transmision  through  the  BS2. 

The  squared  terms  of  the  signal  and  scattered  light  are  of  the  same  form,  i.e., 

£Cw  “ -4C2W  cos2(0)  - j A 2W.  (5) 


Recalling  that  the  scattering  processes  are  assumed  reciprocal  so  that  each  scattered  component 
in  the  CW  direction  is  equal  to  that  in  the  CCW  direction,  then,  by  defining  an  optical  intensity 
/ = \A  |2,  we  get,  for  the  averaged  photodetector  output  current,  the  expression 


i,  ~ D 


1 

2 


lcw  cos(A <f>) 


(6) 


where 


D 


and 


V 

v 

Q 

h 


is  the  detector  conversion  factor  with 
hv 

quantum  efficiency  of  the  detector 
laser  frequency 
electronic  charge 
Planck’s  constant 
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Clearly,  if  ft  = constant,  the  term  that  contains  cos (A<£)  is  a time-independent  signal 
term.  Hence,  the  detection  process  is  necessarily  a dc  mode,  which  can  be  susceptible  to  low- 
frequency  noises. 


Using  a similar  argument,  if  the  similar  detector  system  is  set  at  fringe  pattern  F2,  we 
arrive  at  an  expression  for  the  time-averaged  photocurrent  from  detector  2: 


zr 


«l 


1 


2 (^cw  + ^FP  + ^RF  + ^RB  + ^BO  + ^Bir)  + y I cw  cos  (&<t>) 


(7) 


where 


a,  — power  reflection  coefficient  of  the  BS, 


and 


D - is  the  detector  conversion  factor  for  detector  2. 
hv 


Note  that  the  quantum  efficienty  17 'may  be  different  from  tj  of  detector  1.  The  sign 
change  in  the  last  term  of  Eq.  (7)  is  due  to  complementarity. 


The  total  photocurrent  from  each  detector  is  obtained  by  integrating  the  intensities  of  all 
the  light  components  in  the  fringe  pattern  over  the  small  area  a that  is  seen  by  the  detector. 
Assuming  that  the  optical  waves  are  perfectly  aligned  and  that  the  quantum  efficiency  77  is  uni- 
form across  the  detector  surface,  we  get  total  photocurrents  for  the  detectors  1 and  2 as 


•t\  = Ja  i\da  “ 


— (/cw  + Iff  + /rf  + 7rB  + ^BO  + ^Brr)  ~ y /Cw  COS  (A0) 


— (P Cw  + T’fp  + P RF  + /*RB  + Pm  + T’bJ  ~ y PCw  COS  (A<£) 


(8) 


- D 'l 


— ( p 


(9) 


where  P and  P'  denote  optical  powers  ( P = la).  Here,  P/P'  = aja7  because  / — /'.  If 
a]  — a j,  then  P — P'.  In  (8)  and  (9),  the  first  term  is  the  total  dc  power  and  the  second  is  the 
ac  signal  term.  The  dc  power  will  produce  shot  noise.  Taking  into  account  the  other  possible 
noise  sources,  the  signal-to-noise  ratio,  based  on  the  peak  electrical  power  measured  at  the  load 
resistor  RL  after  the  total  photocurrent  being  amplified  with  gain  G and  filtered  at  baseband  B0 
is  of  the  form  (according  to  standard  homodyne  technique  [11]), 


D'G'P^RJA 


(10) 


2 qG2\ 


2 (^cw  + P rr  + /’rf  + /’rb  + /’bo  + P b„)  + DP b + lD 


+ AkTBn 
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Ph  = externa!  light  induced  background  noise 

/[)  — detector  dark  current  or  other  kinds  of  current  noises  depending  on  the  kind  of 
detector  used 

AkTB„  = the  thermal  noise  term,  where  with  T is  the  absolute  temperature  and  k is 
Boltzmann’s  constant. 

If  the  optical  powers  are  high  and  the  other  noises  can  either  be  low  or  eliminated,  then 
we  are  in  a shot-noise-  or  quantum-noise-limited  situation  whereby  Eq.  (10)  is  simplified  to 

A =1_D_ PA (11) 

N neak  . 4 (^cw  + P FP  + P RF  + P RB  + P BO  + P B„) 


for  detector  1,  and  similarly, 

U - 1 , 7 _ (.2) 

l peak, 2 4 (^cw  + ^FP  + ^*RF  + ^RB  + ^B0  + ^Bir) 

for  detector  2.  Furthermore,  if  the  components  of  the  scattered  light  in  the  optical  fiber  are  all 
neglected,  as  was  done  in  Vali  and  Shorthill  [3,4],  the  peak  S/N  for  a single  detector  reduces  to 

( s)  I DP <=*  I (13) 

Mpeak"  4 *.  “4  hvB„\ 

which  is  one-fourth  of  the  S/N  for  a standard  heterodyne  technique  with  a synchronous  detec- 
tor, and  one-half  of  that  of  the  baseband  direct-detection  technique.* 

Next,  we  discuss  a differential  scheme  for  A<£  or,  more  precisely,  cos  (A <f>)  measurement. 
The  scheme  is  illustrated  in  Fig.  4.  As  we  alluded  to  earlier,  the  advantage  of  this  scheme  is  to 
eliminate  or  reject,  by  subtraction,  the  common-mode  or  completely  correlated  extraneous 
effects  from  the  two  detectors.  In  doing  so,  the  signal  level  will  double  because  of  the 
difference  in  the  signs  of  the  two  ac  terms  in  Eqs.  (8)  and  (9)  above,  while  the  uncorrelated 
noises  will  add  in  the  mean-square  sense.  However,  to  match  or  balance  the  two  detectors  such 
that  the  dc  components  of  the  two  total  photocurrents  will  completely  cancel  in  subtraction,  we 
must  make  the  coefficients  of  the  two  terms  in  Eqs.  (8)  and  (9)  equal.  We  can  accomplish  this 
by  selecting  17,  17',  au  a2,  and  “1  such  that  Dat  — D'a2a t,  yet  still  keep  both  a , and  a2  small  as 
compared  with  the  central  fringe  widths.  After  we  do  that,  the  total  photocurrent  to  be 
amplified  is 

/,  = iT2  — iT\  — Da  | /cw  cos  (A<j>) 

- DP cv  cos  (A</>).  (14) 

Let  the  gain  of  the  differential  amplifier  be  G.  The  peak  signal  power  is,  for  a load  resistor  RL, 

Rl/J-  D2G2PXRl  (15) 

and  the  noise  power,  in  the  quantum-noise  limit,  is 

Pl  'n  “ (//vt  + 'ni^P l^2 

- 2 qB0G2RLD(Pcv+  Fp  + PRf  + PRH  + PB0+  PBn).  (16) 

'See  Ref  It,  Eqs  (10-45)  and  (10-13)  of  Chapter  10. 
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Fig.  4 — Differentia!  detection  scheme  with  balanced  inputs 
to  eliminate  common-mode  noises  adaptable  to  the  optical 
setup  or  Fig.  I 


whence 

(JlI  _ jl  _ _D Pc. 

l ^Ipeak.dilT  <7  cw+*>FP+7,RF  + + (“bJ 


(17) 


Again,  if  scattering  in  fiber  is  neglected,  we  arrive  at 

Ij*|  JL  Pc~ 


L vpc. 

2 hv 


(18) 


Compared  with  Eq.  (13)  for  a single  detector,  the  differential  technique  enables  us  to  gain  a 
factor  of  two  in  peak  electrical  power  S/N  ratio.  Equations  (11)  and  (12)  give  peak  S/N  for  a 
single  detector,  and  Eq.  (17),  for  the  differential  scheme,  all  in  the  quantum-noise-limited 
sense.  They  are  correct  formulas  for  the  optical  configuration  shown  in  Fig.  1 as  used  by  Vali 
and  Shorthill  [4],  indicating  that  their  formulas  would  overestimate  the  sensitivity  by  a factor 
of  4. 


LIGHT  SCATTERING  IN  OPTICAL  FIBERS 


As  seen  from  Eqs.  (11)  or  (12),  and  (17),  quantum-noise-limited  S/N  ratios  depend  on 
those  scattered-light  components  that  are  implicit  functions  of  the  fiber  length  L.  In  order  to 
see  how  they  depend  on  L and  what  their  magnitudes  are,  we  need  to  know,  in  addition  to  local 
scattering  coefficients,  their  distribution  laws  which  would  enable  us  to  calculate  the  cumulative 
powers  at  the  two  ends  of  an  optical  fiber  of  length  L.  In  the  following,  we  shall  show  how  the 
S/N  is  degraded  by  the  added  scattered  light,  based  on  the  best  available  fiber  scattering  loss 
data.  But,  first,  let  us  derive  the  distribution  laws. 

Distribution  Laws 


Direct  Beams 


By  the  direct  beams  we  mean  those  parts  of  the  two  optical  beams  that  have  not  suffered 
scattering  and  absorption  when  propagating  through  the  fiber.  Because  of  the  scattering  and 
absorption,  which  are  the  main  causes  for  the  deletion  of  optical  energy,  the  direct  beam 
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experiences  attenuation.  If  constancy  of  the  scattering  and  absorption  coefficients  and  their  spa- 
tial homogeneity  along  the  fiber  are  assumed,  the  law  of  attenuation  of  the  direct  beam  will  be 
exponential,  i.e., 


where  P,  is  the  power  of  the  beam  at  the  input  end  of  the  fiber,  P(x)  is  the  power  in  the 
transmitted  beam  at  a distance  x from  the  input  end,  and  aT  is  the  total  attenuation  coefficient 
which  is  normally  expressed  as  the  sum  of  the  scattering  coefficient  and  the  absorption 
coefficient.  Let  as  be  the  scattering  coefficient,  and  aabv  the  absorption  coefficient,  then  we 
have  a T — as  + aib5. 

In  applying  Eq.  (19)  to  a realistic  fiber,  we  must  recognize  that  (a)  the  total  attenuation 
coefficient  a r consists  of  contributions  from  intrinsic  scattering  (13  — 15)  and  absorption  as  well 
as  from  extrinsic  scattering  induced  by  external  effects  such  as  bending  and  vibration  of  the 
fiber  (16),  and  that  (b)  the  effect  of  multiple  scattering  (17—19)  has  been  neglected. 

Scattered  Light 

Each  scattered-light  component  can  be  characterized  by  two  parameters;  one  is  the  scatter- 
ing coefficient  a„  where  the  subscript  / indicates  the  Ah  kind  of  scattering  process;  the  other  is 
the  factor  that  quantifies  the  percentage  of  trapping  of  the  scattered  light  in  the  fiber.  We 
assume  that  the  a,’s  are  constant,  pursuant  to  our  above  discussion,  and  that  the  percentage  of 
trapping  remains  unchanged  as  the  scattered  light  propagates  down  the  length  of  the  fiber. 

The  percentage  of  trapping  the  forward-scattered  power  is  denoted  by  F and  that  of  the 
backward-scattered  power,  B.  We  will  derive  the  distribution  laws  for  the  forward  and  the  back- 
ward components  in  the  single-scattering  approximation.  This  approximation  is  justifiable,  in 
that  the  multiple  scattering  at  most  will  result  in  only  a few  percent  of  the  already  small  scatter- 
ing losses. 

Forward-Scattered  Components  — The  symbols  and  coordinates  of  forward  scattering  are 
defined  in  Fig.  5a.  The  unscattered  part  of  the  power  reaching  a station  x from  the  input  end  is 
given  by 

P(x)  - P,e'arI. 

The  power  scattered  from  the  small  element  dx  is 
dP(x ) — — P(x)asdx 

- - P,e~aTXa,dx  (20) 

where  a,  — scattering  coefficient  (—  aT  — «,(*),  and  the  minus  sign  means  a loss.  Then  the 
part  trapped  in  the  forward  direction  is 

dPb r)  - - FP,e~°T 'asdx.  (21) 

When  this  part  reaches  a downstream  station  xh  its  magnitude  is  reduced  by  a factor  e "T>''  r> 
So, 

dP,(x)  - - FP,e  " T*e  “ r<Jt|  _Jr,a  ,<fx 

- - FP,ay°TX'dx. 
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If  we  integrate  dP5  from  0 to  x,,  we  obtain  all  the  scattered  light  produced  by  the  fiber  from  the 
input  end  up  to  the  station  xj.  The  integrated  result  is 

/>,(*,)  - - FP,a,e~°'TXlxl.  (22) 

When  xx  - L,  i.e.,  the  output  end  of  the  fiber,  we  then  have  the  total  forward-scattered  light 

P,(L)  - - FP,a,e~aTLL.  (23) 

This  is  the  distribution  law  for  the  forward-scattered  light. 


Backward-Scattered  Components  — According  to  Fig.  5b,  the  scattered  power  from  the 
small  element  dx  back  to  a station  x,  with  B fraction  trapped  is 

dP,(x)  - - BP;e~aTXa,e'arix~Xl)dx 

- - BP,a ,earXle~la rXdx.  (24) 

Then  the  total  backscattered  light  reaching  X)  from  the  length  of  fiber  X)  — L is 
/>,(x.)  - - BPia,earX'fxL  e-2aTXdx 


_ BP)  a,  ^-aTOL-xt)  _ ^-orx,j 

2 «r 


(25) 


Hence,  the  total  backscattered  light  arriving  at  the  input  end  is  given  by  letting  x,  — L in 
Eq.  (25): 


P,(0) 


.SLl 


2 aT 

This  is  the  distribution  law  for  the  backscattered  light. 


(26) 
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Fig.  5 - Definition  of  coordinates  and  scattering  fiber  ele- 
ments of  the  forward  scattering  (a)  and  backward  scattering 
(b)  processes  for  the  derivation  of  scattered  light  distribu- 
tion laws. 


12 


T 


NRL  REPORT  8250 


l 


i 

. 

j 

! 


Trapping  Factors 

For  the  fundamental  mode  (HEM)  in  a single-mode  fiber,  the  power  is  distributed  into  the 
cladding  [20,211.  Therefore,  when  we  discuss  about  the  scattered  light,  in  general,  we  should 
divide  the  contribution  into  two  parts:  one  from  the  inhomogeneities  located  randomly  inside 
the  core,  and  the  other  from  those  in  the  cladding.  Rays  scattered  from  inside  the  core  that 
strike  the  core— cladding  interface  at  an  angle  larger  than  the  critical  angle  defined  by  9cr  — 
sin ~l(nCi,d/« core)  will  be  trapped  inside  the  core.  Rays  striking  at  an  incidence  angle  0 , < 0cr 
will  escape  into  the  cladding  and  will  be  either  trapped  or  further  escape  into  the  medium  sur- 
rounding the  cladding,  depending  on  whether  the  index  of  refraction  of  the  surrounding 
medium  is  smaller  or  larger  than  that  of  the  cladding,  respectively.  Some  rays  initiating  from 
scattering  centers  in  the  cladding  may  be  trapped  in  the  cladding  if  its  index  of  refraction  is 
larger  than  that  of  the  surrounding,  whereas  others  may  refract  into  the  core  and  be  trapped 
into  a helical  course.  Complications  arise  if  we  take  into  account  real  effects  of  the  intensity 
distribution  of  the  HE(I  mode,  mode  distortion  due  to  bends  in  the  fiber,  and  the  core- 
cladding  interface  irregularities.  We  will  limit  our  discussion  to  a simple  case  where  we  con- 
sider only  the  more  significant  scattering  components  arising  from  the  inhomogeneities  distri- 
buted in  the  core  and  the  core— cladding  interface.  The  computation  of  trapping  factors  for 
these  components  will  follow  a simple  approach  by  J.  Stone  [22],*  who  neglected  the  scattered 
light  contributed  by  the  cladding.  This  could  be  justified  by  the  fact  that  the  intensity  of  the 
direct  beam  in  the  cladding  is  much  lower  than  that  in  the  core  [20].  To  demonstrate  the  possi- 
ble influence  of  all  the  trapped  scattered  light  in  the  fiber,  we  will  consider  the  case  where  no 
mode-stripping  or  index-matching  compound  is  applied  on  the  cladding,  together  with  the  more 
favorable  case  where  mode  strippers  are  used  to  reduce  the  cladding  trapping. 


Since  the  percentage  of  the  trapped  light  determines  the  magnitude  of  the  influence  of  the 
scattered  light  on  the  ultimate  sensitivity  of  the  detection  system,  we  will  assess  the  magnitudes 
of  the  trapping  factors,  F and  B.  We  must  note  that  if  there  is  no  trapping,  or  a very  small 
amount  of  trapping,  then  the  contribution  of  scattered  light  comes  only  from  both  ends  of  the 
fiber  instead  of  from  an  integration  of  all  the  trapped  scattered  light  along  the  whole  length  of 
the  fiber. 


Rayleigh  and  Brillouin  Scattering 


According  to  J.  Stone  [22],  the  formula  for  computing  the  one-way  fraction  of  trapped 
power,  nlrap,  for  a bare  cladding  in  air  is 


ntrip  - f0  d<t>  fQ  ° (1  + cos2  9)  sin  9 d9 
_ 4ir  4 |_1 ill 

2 3 J ncore  3ec3oreJ 


(27a) 


*A  more  refined  version  of  ihe  refraction  and  reflection  at  the  core— cladding  interface  was  given  more  recently  by  I P 
Dakin  and  W A Gambling,  Opt.  Commun  10.  1 95- 1 98  (1974).  Because  we  arc  interested  in  the  firsi-order  effects,  we 
have  chosen  Stone's  formula  for  symmetric  scattering  Dakin  and  Gambling's  formula,  which  lakes  into  account  the 
effects  of  depolarization,  transverse  intensity  distribution  in  Ihe  core,  and  the  finite  angular  width  of  the  propagating  in- 
terface  will,  of  course,  be  used  when  more  refined  calculations  are  needed 
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where  #n  cos  1 ( 1 / /»tore)  is  ihe  limiting  scattering  angle  below  which  a ray  will  be  trapped  in 
the  cladding,  /»UIIC  is  the  index  of  refraction  of  the  core  material,  and  the  factor  (1  + cos2#)  is 
the  angular  factor  for  Rayleigh  scattering  for  unpolarized  light.  The  first  term  in  Eq.  (27a)  is 
the  fraction  of  the  total  Rayleigh  scattered  power.  It lo,.|(  = 16ir/3. 


For  a cladding  with  mode  stripper  which  has  an  index  of  refraction  higher  than  that  or  the 
cladding,  then,  light  rays  with  scattering  angle  close  to  #„  will  escape.  The  amount  that  will  be 
trapped  is  determined  by  the  internal  refraction  at  the  core— cladding  interface,  i.e. 
#„  - cos  nrc),  and  Equation  (27)  becomes 


2 w ntl 

I,  Jo  * 1 + cos2#)  sin  # (/# 


— — 

4 

22  clad  "clad 

2 

3 

22  wire  3 nj„rc 

(27b) 


Therefore,  the  trapping  factor  F is 


for  bare  cladding,  and 


f = ^ Irap  ^ 
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(28a) 


(28b) 


for  a mode-stripped  or  index-matched  cladding.  Because  of  symmetry  of  the  Rayleigh  scatter- 
ing, the  backward  trapping  factor  B is  equal  to  F.  According  to  J.  Schroeder,  el.  al.  (23],  Bril- 
louin  scattering  has  the  same  spatial  angular  distribution  as  the  Rayleigh  scattering.  The  trap- 
ping factors  for  Brillouin  scattering  are,  therefore,  identical  to  the  Rayleigh.  For  given  materi- 
als of  core  and  cladding,  the  trapping  factors  can  be  evaluated.  For  example,  let  us  consider  an 
optical  fiber  with  fused  silica  core  (#c„rc  - 1.4585)  and  a borosilicate  glass  cladding 
foSiOjilBjO,)  with  an  index  of  refraction  about  0.3%  less  than  that  or  pure  fused  silica  (at 
X - 1.06  Aim,  neM  - 1.4541).  The  limiting  scattering  angle  #,,  = 46.7°  and  4.44°  for  the  bare 
cladding  and  index  matched  cladding,  respectively.  Based  on  Eq.  (28a),  F = B - 0.2  and  2.3  x 
10  ’,  corresponding  to  the  respective  #„.  We  can  immediately  draw  the  conclusion  that  the 
fiber  to  be  used  in  a ring  interferometer  should  be  equipped  with  mode  strippers,  or  index- 
matching components,  to  reduce  the  trapping  factors  by  about  100  times  for  the  Rayleigh  as 
well  as  the  Brillouin-scattered  components. 


Mie  and  Forward  Peak  Seauering 

In  the  study  of  light  scattering  in  the  optical  fibers,  a strong  forward  scattering  component 
has  been  observed  both  for  a single  mode  fiber  (24]  and  a multimode  fiber  (251.  These  peaks 
were  observed  just  outside  the  cladding  when  the  fibers  were  immersed  in  index-matching  fluid. 
This  peak  cannot  be  explained  by  Rayleigh  scattering  theory  and  is  generally  thought  to  be 
caused  by  two  factors:  namely,  the  Mie  scattering  and  the  scattering  due  to  core— cladding 
interface  corrugation.  Mie  scattering  was  first  studied  by  Rawson  [26]  for  a bulk  glass  with 
imbedded  spherical  particles.  These  micron-size  particles  are  the  homogeneities  frozen  in  dur- 
ing solidification  (27).  However,  inasmuch  as  these  spherical  particles  are  liable  to  be  drawn 
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into  filaments  aligned  with  the  fiber  during  the  fiber-drawing  process,  Rawson  devised  a theory 
for  light  scattering  due  to  randomly  distributed,  but  aligned,  filaments  in  a fiber  (28] . On  the 
other  hand,  the  interface  scattering  process  was  first  studied  theoretically  by  Marcuse  |29). 
Since  then,  there  have  been  more  studies  [30—33],  but  none  is  experimental,  presumably 
because  of  the  difficulties  involved  in  measurements.  Experimentally,  it  is  very  difficult  to 
separate  these  two  contributions  [24,25].  Nevertheless,  experimental  results  do  establish  the 
ratios  of  the  total  forward-scattering  component  an  the  Rayleigh  component  after  other  leaky 
modes  are  stripped.  We  will  use  these  ratios  to  calculate  the  forward-peak  scattering  com- 
ponent from  the  knowledge  of  the  Rayleigh  scattering 

It  is  also  conceivable  that  these  losses  are  dependent  on  the  fiber  drawing  processes, 
hence  are  variable  from  case  to  case.  But,  in  order  to  render  some  degree  of  generality  for 
analysis  sake,  we  will  assume  that  the  ratios  mentioned  above  will  stay  constant  from  fiber  to 
fiber,  provided  that  the  fibers  are  drawn  by  a similar  process  from  similar  preforms. 

Our  current  knowledge  of  a fiber's  ability  to  trap  the  forward-peak  scattering  is  too  limited 
to  enable  us  to  quantitatively  assess  the  trapping  factor  for  a general  case.  For  example,  the 
question  of  how  much  of  the  forward  peak  will  be  guided  by  the  fiber  when  heavily  mode 
stripped  does  not  have  an  answer,  so  it  is  with  the  question  about  the  effects  of  fiber  bending 
on  the  propagation  of  the  peak,  for  all  the  measurements  so  far  have  been  for  a short,  straight 
section.  However,  because  a strong  forward  peak  is  observed  outside  of  the  fiber  when 
immersed  in  an  index-matching  fluid,  we  see  clearly  that  proper  mode  stripping  would  eliminate 
at  least  a large  portion  of  the  forward-peak  scattering  component.  To  demonstrate  how  high 
the  forward  peak  would  be  if  fully  trapped,  we  will  discuss  this  case  in  detail.  This  will  be  fol- 
lowed with  a discussion  of  a case  wherein  all  the  forward  peak  is  assumed  to  have  been 
removed  from  the  fiber. 

The  trapping  factors  for  the  former  case  are  F = 1 and  B = 0,  whereas  the  trapping  fac- 
tors for  the  latter  case  are  F — 0 and  B = 0.  Of  course,  in  the  latter  case,  we  no  longer  speak 
of  a distribution  law.  Instead,  we  should  consider  the  part  that  is  scattered  from  the  exit  face  of 
the  fiber,  which  depends  on  the  end-surface  conditions. 

Scattering  Coefficients 

As  a basis  for  our  discussion,  we  choose  a fiber  that  has  a total  loss  rate  of  4 dB/km  at 
X = 0.633  /im.  We  further  assume  that  the  total  loss  rate  is  due  purely  to  scattering,  in  that 
the  absorbing  centers  in  the  glass  fiber  [34]  can  be  completely  removed  with  advanced  tech- 
niques. In  other  words,  we  are  considering  here  a fiber  that  has  a loss  mechanism  in  the  fre- 
quency range  of  interest  due  to  scattering  alone.  We  will  use  a weakly  guiding  fiber  as  defined 
by  Gloge  [21]  that  has  a fused  silica  core  for  quantitative  assessment  of  scattering  coefficients. 

Brilloum  Scatterinn  Coefficient 

Rayleigh  and  Brillouin  scattering  in  glass  have  been  studied  extensively  in  the  past 
[23,35  — 38],  But,  it  seems  that  no  detailed  study  of  Brillouin  scattering  from  an  optical  fiber  has 
been  reported  in  the  literature,  except  the  work  by  Rich  and  Pinnow  [39],  which  reported  the 
measurements  of  Brillouin  scattering  from  both  the  core  and  the  cladding  of  multimode 
borosilicate-pure  fused  silica  waveguide,  and  that  of  E M.  Dianov,  el  al.  [40] . Brillouin  scatter- 
ing from  a single-mode  fiber  is  still  alcking.  In  view  of  this,  we  will  use  formulas  for  bulk 
material  to  calculate  Brillouin  scattering  for  the  single-mode  fiber 
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According  to  Rich  and  Pinnow,  the  Brillouin  scattering  coefficient  is  given  by  the  Formula, 


where  A is  the  optical  wavelength  in  Free  space,  k is  Boltzmann’s  constant,  T is  the  absolute 
temperature,  n is  the  index  oF  reFraction  oF  the  glass,  p!2  is  the  applicable  photoelastic  tensor 
component,  p is  the  density,  and  F'is  the  longitudinal  mode  acoustic  velocity.  For  Fused  silica, 
they  give  0.40  dB/km*  For  aSB  at  room  temperature  and  A -0.5145  pm.  According  to  J. 
Schroeder,  et  al.  [23),  the  Brillouin  scattering  coefficient  can  be  obtained  through  the  Rayleigh 
ratio  evaluated  at  90°,  i.e.,  aSB  - ( &ir/3)R90°  For  linearly  polarized  light,  where  the  Rayleigh 
ratio  Aw* is 

^ (/»Vi22)  kT  (pk'2)-'.  (30) 

Equation  (30)  gives  an  aSB  identical  to  Eq.  (29).  IF  we  use  the  thermophysical  and  optical  data 
oF  Si02  provided  by  J.  Schroeder,  et  al.  [23]  in  Tables  1—3  namely,  V — 5.92  x 105  cm/s, 
P- 2.211  g/cm3,  p,2  - 0.286,  T - 293K,  n — 1.458,  and  A - 0.633  pm,  we  get  #90-  — 
0.538  x 10'7  cm1;  hence  aSB  — 0.195  dB/km*  For  the  V—  V polarization  component.  For  the 
unpolarized  light  under  consideration,  aSB  is  twice  as  large:  aSB  — 0.39  dB/km,  which  compares 
well  with  the  value  0.4  dB/km  given  by  Rich  and  Pinnow  For  A - 0.5145  pm.  In  what  Follows, 
we  will  use  the  value  0.40  dB/km  as  the  intrinsic  Brillouin  scattering  loss  For  the  fiber  at 
A - 0.633  pm. 

Effective  Forward-Peak  Scattering  Coefficient 

From  experimental  results,  Rawson  [24]  established  that.  For  the  particular  single-mode 
fiber,  out  of  the  total  scattering  loss  rate  oF  about  10  dB/km  at  A - 0.633  pm,  the  Forward  peak 
took  a share  oF  2.5  dB/km,  while  the  Rayleigh  scattering,  the  balance,  about  7.5  dB/km.  The 
ratio  oF  the  Forward  peak  over  the  Rayleigh  scattered  light  is  about  1:3.  In  the  same  paper, 
Rawson  also  established  that  the  spectral  dependence  oF  the  Rayleigh  component  Follows  the 
1 /a4  law  within  7%.  The  results  measured  at  A - 1.06  pm  show  a total  loss  oF  about  1.05 
dB/km,  with  Forward  peak  about  0.17  dB/km  and  Rayleigh  component  about  0.88  dB/km.  The 
ratio  is  about  1:5,  instead  oF  1:3.  It  is  interesting  to  note  that  the  Forward  peak  may  not  obey 
the  1/A4  law,  indicating  that  the  Forward  peak  could  be  caused  by  effects  other  than  the  Ray- 
leigh scattering. 

The  same  kind  oF  measurement  was  again  careFully  done  by  Reeve,  et  al.  [25)  with  a mul- 
timode fiber.  They  reported  a residual  Forward  peak  constituting  about  one-fiFth  oF  the  Ray- 
leigh scattered  power  at  0.633  pm,  and  attributed  the  cause  to  large-scale  (w.r.l.  the 
wavelength)  imperfections.  They  also  Found  the  1/A4  characteristics  associated  with  the  scat- 
tered loss  aFter  the  Forward  peak  had  been  subtracted  out.  ThereFore,  in  our  calculations,  we 
will  adopt  a value  oF  1:5  For  the  ratio  of  the  forward  peak  power  over  that  of  the  total  Rayleigh 

’Rich  and  Pinnow  did  not  give  the  polarization  slide  of  the  light  when  giving  Inc  value  0 4 dB/km  Our  calculations 
based  on  Ihcir  formula  dig  I)  and  Schrocdcr's  result  in  values  only  one-hall  of  it  il'  linear  polarization  is  assumed  We 
concludc  that  0 4 dB/km  would  be  correct  lor  unpolari/cd  light,  which  is  what  we  presume  now 
tl  cm'1  — 4.33  x I05  dB/km  is  the  conversion  factor 
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component.  Using  this  ratio,  we  are  able  to  ascribe  an  effective  scattering  coefficient  for  the 
forward  peak  from  the  Rayleigh  scattering  coefficient. 

Partition  of  Scattering  Losses 

Since  Brillouin  scattering  is  intrinsic,  we  will  let  it  be  constant.  The  total  loss  is  assumed 
to  be  4 dB/km;  the  Rayleigh  scattering  and  the  forward  peak  scattering  are,  then,  partitioned 
according  to  the  formulas 

“kp  “ (1/5)  aSn 
and 

«sk  “ «/■  “ “SB  - “i  p (31) 

where  u(P  is  the  equivalent  scattering  coefficient  for  the  forward  peak,  which  is  one-fifth  of  the 
total  Rayleigh  scattering  components  as  we  established  above.  Since  a , ■*  4 dB/km  = 9.24  x 
10  7 cm  1 and  aSH  - 0.4  dB/km  - 92.4  x 10  7 cm  1 for  X = 0.633  (im,  therefore,  by  Eq. 
(31),  we  get,  t»j  P — 13.86  x 10  7 cm  1 and  «Sr  = 69.3  x 10  7 cm  '.  Also,  «Sr/“  / = 0.1  and 
uIP/or  / *■  0. 1 5. 

For  comparison,  we  also  calculate  these  scattering  coefficients  for  X = 1.1  pm.  The  rea- 
sons for  choosing  this  wavelength  will  become  clear  later  on.  We  use  the  fact  that  the  Rayleigh 
scattering  coefficient  scales  as  1/X4,  and  the  fact  that  the  forward  peak  scales  as  1/X3,  as  esta- 
blished by  Rawsom  124,281.  Also,  according  to  Rawson’s  measurements  (24]  for  = 1.06  pm, 
the  forward  peak  is  about  0.17  dB/km,  whereas  the  Rayleigh  is  0.88  dB/km.  The  ratio  ol  Ray- 
leigh scattering  over  the  forward  peak  is  5.18.  So,  for  «•  1.1  this  ratio  is  about  5.0. 
Therefore,  the  forward  peak  is  still  one-fifth  of  Rayleigh  component.  The  Brillouin  scattering 
coefficient  also  scales  as  1/X4,  as  can  be  seen  from  Eq.  (29).  Its  value  for  X = 1.1  /^m  is  now 
0.043  dB/km.  To  calculate  all  the  other  scattering  components,  we  assume  that  the  ratio 
between  the  Brillouin  and  the  Rayleigh  components  remains  unchanged  when  the  wavelength  is 
changed  from  0.633  pm  to  1.1  pm\  i.e.,  <xSH/asl(  - 1/7.5.  Hence,  in  conjunction  with  the  con- 
dition «//.  — (1/5)«SK,  we  get  a r ~ 0.431  dB/km  and  «KP  — 0.065  dB/km  for  «SK  — 7.5  x 
uSB  — 0.323  dB/km.  The  percent  ratios  are:  uSH/«/  “01.  «sr/“/  ^ 0.751,  and  aPP/a , — 
0.149.  The  value  of  a,  represents  a realistic  assessment  of  the  intrinsic  scattering  loss.  It  is 
approached  by  the  current,  rapidly  advancing  fiber  manufacturing  technology.2  We  also  note 
that  the  total  loss  rate  does  not  scale  as  X 4 if  the  cladding  mode  is  not  stripped2  and  if  the 
forward  peak  contributes  significantly. 

Because  the  index  of  refraction  of  the  fused  silica  changes  very  little  in  the  wavelength 
range  0.5  - 1.3  riti  [41,42],  the  fiber  trapping  factors,  as  a weak  function  of  the  index  of 
refraction  of  the  core,  nL„,c,  do  not  change  appreciably.  Therefore,  we  will  use  the  same  trap- 
ping factors  for  both  wavelengths. 

For  readers’  convenience,  we  conclude  this  chapter  by  summarizing  the  useful  informa- 
tion in  Table  I. 
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Table  1 — Summary  of  Light  Scattering  Information  in  a Typical  Fused  Silica  Fiber 
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RELATIVE  POWER  DISTRIBUTION 

AND  SIGNAL-TO-NOISE  DEGRADATION 

With  ifite  scattering  losses  and  trapping  factors  determined  above,  we  can  now  evaluate 
numerically  the  relative  power  distribution  laws  of  all  the  light  components  in  a fused  silica 
fiber  of  length  L.  Also  evaluated  will  be  the  degradation  factor  of  the  S/N  ratio  associated  with 
fiber  scattering  losses. 


Relative  Power  Distributions 

The  relative  power  distribution  laws,  for  fiber  with  bare  cladding,  lake  the  following  forms 
(one  way  only): 

PJU  = PCJ\\  - e "'1 


P Kp/d 


1 

J“SK 


Le 


r> 


/V/n  - '0.2 a^)Le"r 


PRH/l\  = (0.1aSK/ar)  (1  - e '2"r^ 


PvJ II  - (0.2aSB)  Le  "r 


PHJ II  = (0.1«SB/ar)ll  - e~2"r'l,  (32) 

where  L is  the  total  length  of  the  fiber  coil,  II  is  a normalization  factor,  II  = ( 1/4) ( I - a^PoC. 
As  described  in  Chapter  1,  the  decimal  numerical  coefficients  are  trapping  factors.  The  reader 
is  reminded  of  the  fact  that,  in  Eq.  (32),  the  backscaltered  components  PKB  and  Pu„  are  due  to 
the  wave  traveling  in  the  opposite  direction  to  that  of  the  other  components.  The  values  of 
these  scattering  coefficients  are  given  in  column  4,  Table  1,  for  both  X =0.633  /am and  A — 1.1 
Mm.  We  have  evaluated  Eq.  (32)  as  a function  of  L.  The  results  are  plotted  in  Fig.  6.  We  see 
that  the  forward-scattering  components  first  increase  with  L and  then  reach  maxima  at  L = 
Lmt%.  After  that  they  decay  continuously,  whereas  the  backscattered  components  increase  con- 
tinuously with  L and  tend  to  saturate  for  large  L.  The  direct  signal  beams  are,  of  course, 
straight  lines  in  the  semilog  plot.  At  L » LmdK,  the  background  noise  will  be  contributed 
mainly  by  backscattered  light.  The  Brillouin  components  are  a factor  of  7.5  below  the  Rayleigh 
and  forward  peak  components  (the  forward  peaks  coincide  incidentally  with  the  Rayleigh  for- 
ward components).  Since  the  total  scattering  loss  rate  for  X = 1.1  /urn  is  about  10  times  less 
than  that  for  X = 0.633  /im,  the  direct  beam  is  obviously  decaying  much  more  slowly.  The 
scattered  components,  in  this  case,  change  with  L at  a much  slower  rate,  although  interestingly, 
they  saturate  to  the  same  levels  at  large  L as  they  would  for  X = 0.633  /zm.  We  see  also  that 
by  going  to  a longer  wavelength,  we  gain  a wider  range  of  fiber  length  within  which  we  can 
achieve  certain  specific  sensitivity  requirements  by  using  10-  or  20-km  lengths  of  fiber  without 
large  degradation  in  S/N  ratios.  We  conclude  that,  from  the  standpoint  of  scattering  losses,  we 
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Eig.  6 — Distribution  laws  for  all  the  scattering  components  considered  for  a 
liber  with  attenuation  coefficients  <r r - 0 43  dB/km  at  A - I I /ini  (solid  lines) 
and  it r = 4 dB/km  at  A = 0 633  fim  (broken  lines).  Due  to  attenuation,  the 
forward  scattered  light  components  rise  first  as  liber  length  is  increased,  reach 
maxima,  and  then  decay,  whereas  the  backscattercd  components  rise  continu- 
ously to  saturation  The  direct  beam  attenuates  much  more  slowly  in  the  case 
A = I 1 /cm  than  that  for  A = 0.633  /iin,  while  backscattercd  noises  saturate  to 
about  the  same  level. 

should  prefer  1.1  /am  to  0.633  |tm  as  the  light  source  wavelength.  With  mode  stripping,  all  the 
scattered-light  components  will  have  magnitudes  lower  by  about  a hundredfold  except  the  for- 
ward peak,  whose  reduction  due  to  mode  stripping  is  not  absolutely  known. 

Signal-to-Noise  Degradation 

To  evaluate  the  degradation  of  S/N,  we  take,  as  an  example,  the  single-detector  case. 
The  expression  for  S/N  (from  Eq.  (ID)  is 

_S  = l _g PL 

N peak  4 hvB0  P c„  + PFf  + P RF  + Prb  + P B0  + Pa* 

- Kf(L)  (33) 

where  K s (\/4)r)/hvB0  is  a constant  independent  of  the  fiber  length,  £,  and  f(L)  — 
P^Ji  1 + P\/Pcw)  <s  a complicated  function  of  Z.,  with  PN  the  sum  of  the  five  scattered  com- 
ponents. 

Physically,  f(L)  represents  the  S/N  ratio  degradation  factor  due  to  light  scattering  in  the 
optical  fiber  with  length  L , and  KP cw  is  the  S/N  ratio  if  scattered  light  components  are 
neglected  as  did  in  Vali  and  Shorthill  (3).  With  the  distribution  laws  as  given  by  Eq.  (32) 
above,  we  have  evaluated  f(L ) = (S/Nlp^/K.  The  results  are  plotted  in  Fig.  7.  Curve  1 is 
for  A = 1.1  /am  and  Curve  2 is  for  A — 0.633  ^m.  The  degradation  of  S/N  in  the  case  A = 
0.633  (im  can  by  no  means  be  neglected.  The  degradation  in  the  case  A = 1.1  jem  is,  by  con- 
trast, very  mild.  This  is  another  reason  for  adopting  a longer  wavelength.  With  mode  strip- 
ping, of  course,  the  degradation  is  much  less  serious. 
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L,  km 


I T(4  7 — Signal-lo-noise  degradation  due  ■<>  til- 
lered light  components  in  a liber  us  a I'unclinn  of 
optical  liber  length  I.  For  a perfect  liber  with  /cro 
loss,  the  S/N  is  purely  signal  quantum-noise  limited 
and  is  equal  to  I here.  Curve  (I)  is  for 

K » 1.1  nm  (043  dB/km):  Curve  (2).  for 
K -0633  jam  (4  dB/km)  So,  at  longer 

wavelength,  the  degradation  is  much  less 
significant 
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PROBLEM  AREAS  AND  THEIR  PROMISING  SOLUTIONS 

From  the  preceding  analyses  and  discussions,  three  problem  areas  have  emerged,  mainly 
due  to  the  stationary  nature  of  the  proposed  optical  arrangement  and  the  homodyne  detection 
schemes. 

Need  for  Stable  DC  Phase  Bias 

The  first  problem  area  is  associated  with  the  fact  that  the  ac  term  of  the  photodetector 
current  (i.e.  the  signal)  is  proportional  to  cos  (A<£);  e g.,  see  Eq.  (8).  This  function  has  least 
sensitivity  at  small  values  of  A0,  since  the  sensitivity  is  proportional  to  t H.Jd(b4> ).  As  seen  in 
Fig.  8,  the  central  fringe  intensity  has  a maximum  sensitivity  in  phase  shift  at  A<£  = rr/2  or 


Fig  8 — The  variation  ol  the  photodciccior  current 
as  a function  of  the  phase  difference  between  the 
two  signal  beams  Maximum  sensitivity  points 
correspond  to  the  points  on  the  photodclccior 
current  curve  where  its  slope  has  maximum  values 
Therefore,  lor  small  ± A»/».  we  need  a dc  phase 
bias  ol  rr/2,  in  order  to  ensure  maximum  sensitivi- 
ty 
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3 7r/2.  In  order  to  ensure  the  location  of  A <t>  at  these  points,  we  must,  therefore,  provide  a con- 
i stant  dc  phase  bias.  Sinee  there  are  two  oppositely  travelling  waves,  we  must  use  nonreciprocal 

elements  based  on  the  well-known  Faraday  effect  143 ] to  introduce  phase  biases.  The  Faraday 
effect  has  been  well  demonstrated  in  magneto-optical  components  such  as  optical  isolators. 

Very  compact  ones  have  been  reported  either  for  the  visible  (44 J or  infrared  (451  region  of  the 
optical  spectrum.  While  many  magneto-optical  materials  (46|  are  available,  the  most  promising 
one  is  represented  by  yttrium  iron  garnet  (YIG,Y3Fe5Ou)  which  has  good  transparency  in  the  s 

l-4-/xm  region  and  high  Faraday  rotation  (=  250  deg/cm  @1.1  gm).  Isolators  without  polar- 
izers (45°  Faraday  rotator)  constitute  the  basic  elements  that  provide  the  90°  dc  phase  bias  for 
the  ring  interferometers.  We  can  obtain  a 45°  Faraday  rotation  by  using  an  approximately  2.3- 
mm-long  YIG  crystal  at  1.25  /xm  wavelength.  Even  smaller  absorption  loss  and  anomalously 
larger  Faraday  rotation  than  the  YIG  have  been  measured  by  Takeuchi  et  al.  (47 1 in  a single 
crystal  of  bismuth-substituted  gadolinium  iron  garnet  (Gd2  ,Bi^Fe30(2  ((  = 0-1.4)).  For 
example,  they  report  Faraday  rotation  of  more  than  3 x 10!  deg/cm  at  X = 1.1  /xm  in 
Gd2Gi|Fe50i2.  with  even  less  absorption  loss  than  YIG.  Such  large  specific  rotation  makes  it 
very  promising  to  develop  miniaturized  components  for  compact  interferometers.  The  availa- 
| bility  of  such  good  characteristics  in  these  garnets,  coupled  with  fiber  loss  as  low  as  0.5  dB/km 

at  1.1  /xm,  the  newly  developed  GalnAsP/InP  double  heterostructure  diode  lasers  emitting  at 
1.1  /xm  at  room  temperature,  reliable  cw  operation  [48],  and  Ge  avalanche  photodiodes  with 
50%  quantum  efficiency  for  X = 0.6-1.65  /xm  [12],  indicate  clearly  the  desirability  of  using  long 
wavelengths  in  ring  interferometric  gyros.  Similar  considerations  have  prompted  Kimura  and  i 

Daikoku  149]  to  propose  a low-loss  1.0-  to  1.4-/xm  region  for  future  optical  fiber  communica- 
tions systems. 

Vulnerability  to  Low-Frequency  Noises 


The  second  problem  is  due  to  the  dc  nature  of  typical  optical  detection  scheme  discussed 
above.  As  such,  it  is  vulnerable  to  those  very  low-frequency  noises  associated  with  the  detec- 
tors and  electronics,  because  we  cannot  subtract  them  out  through  the  differential  schemes  due 
to  their  lack  of  correlations.  These  low-frequency  noises  include  I //-type  noise,  generation  and 
recombination  noise,  partition  noise,  etc.,  with  1 / /-type  noise  the  most  detrimental  [50,51], 
Therefore,  we  will  only  discuss  the  1// noise  effect. 


The  1// noise  is  common  to  all  solid-state  devices.  It  is  a type  of  noise  that  has  a 1 If' 
power  spectrum,  where  0.8  < n < 2.0,  with  « ~ 1 the  most  probable  value.  The  1//  noise 
spectrum  has  extremely  low  frequencies;  the  high-end  cutoff,  however,  can  vary  from  less  than 
1 Hz  to  tens  of  kHz,  depending  on  the  detectors.  To  illustrate  this,  we  show  in  Fig.  9,  two  typ- 
ical spectra  of  1//  noises,  the  upper  one  shows  I//  noise  associated  with  the  InAs  photodiode, 
as  measured  by  J.  Hanlon  and  S.F.  Jacobs  (52],  and  the  lower  one  that  of  the  RCA  1P28  pho- 
tomultiplier, as  obtained  by  C.  Smit  el  al.  (53].  We  see  that  (a)  the  photomultiplier  and  the 


InAs  photodiode  both  have  wideband  I//  noises,  (b)  at  low-frequency  ends,  they  are  orders  of 
magnitude  higher  than  the  shot  noise,  and  (c)  the  photomultiplier  shows  a very  low  cut-off  fre- 
quency at  the  high  end  (<  0.1  Hz)  of  (he  spectrum.  It  may  be  noted  that  al  / ~ 0.3  Hz,  the 
magnitude  of  the  I//  noise  power  is  of  the  same  order  as  that  of  the  shot  noise. 

We  may  conclude  that  the  1//  noise  is  not  too  serious  if  the  photomultiplier  is  used  in 
lieu  of  other  solid  state  detectors  in  dc  or  homodync  detection  schemes,  provided  that  the  size 
of  the  photomultiplier  can  be  accommodated  in  a system  and  that  the  measurements  do  not 
involve  an  average  time  longer  than  20  s. 
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Kig.  9 — Typical  solid-slalc  dclcclor  noise  spectrum  (upper 
figure.  InAs,  Ref.  52).  and  photomultiplier  noise  spectrum 
(lower  figure,  Ref.  53)  showing  the  effect  of  I//  noise  at 
low  frequency  ends. 

However,  if  we  demand  that  the  size  of  a ring  interferometer  rate  gyro  be  as  compact  as  a 
practical  ring-laser  gyro,  then  we  must  select  an  avalanche  photodiode  to  replace  a photomulti- 
plier because  an  avalanche  photodiode  can  be  much  smaller  in  physical  size  and  has  the  same 
internal  high  gain  as  the  photomultiplier.  In  general,  the  nonstationarity  of  the  statistics  of  the 
1//  noise  (54]  does  not  permit  the  long  averaging  time  often  required  in  slowly  varying 
phenomena.  Under  these  circumstances,  the  discussed  optical  scheme  will  fail  to  give  a high 
degree  of  sensitivity.  This  leads  to  the  third  problem. 


Need  for  High  Sensitivity 

The  third  problem  is  the  need  to  achieve  a high  degree  of  sensitivity.  To  this  end,  we 
must  modify  the  optical  scheme  either  to  incorporate  some  kind  of  modulation  technique  into 
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the  homodyne  setup  or  to  adopt  a heterodyne  setup.  The  purpose  ol"  both  modifications  is  to 
defeat  the  I//  noise  by  operating  the  system  at  a higher  frequency  (carrier  frequency  is  now 
shifted  away  from  the  dc  limit)  to  ensure  shot-noise-limiled  performance. 


In  the  case  of  heterodyning,  Hanlon  and  Jacobs  152)  have  demonstrated  the  feasibility  of 
using  a bandwidth  as  narrow  as  0.5  Hz  at  A = 1.15  gm  with  a 50-kHz  phase  modulation  to 
achieve  the  shot-noise-limiled  S/N  ratio.  A true  heterodyne  scheme  also  has  the  capability  of 
providing  a strong  local  oscillator  beam  to  further  improve  the  S/N  ratio,  f igure  10  shows  one 
of  the  heterodyne  schemes  that  we  consider  applicable.  It  has  a local  oscillator  beam  (fre- 
quency) which,  alter  passing  through  the  first  acoustic  modulator  (A  M ),  where  it  is  modulated 
at  frequency  to2,  is  split  into  a zeroth  order  beam  and  a first-order  beam  at  ai(  - a»?,  being 
the  laser  frequency.  The  unshifted  zeroth  order  beam  proceeds  to  the  second  acoustic  modula- 
tor, modulated  at  to,,  to  further  generate  a Bragg-diffracted  beam  at  w,  — <u,.  Optical  mixing 
occurs  at  the  surfaces  of  detectors  A),  and  D7.  The  demodulated  signals  are  then  sent  to  a 
phase  meter  to  produce  a direct  phase  difference  reading  of  2A<|>. 
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t ig  10  — An  example  of  ac  detection  scheme  lor  improve- 
ment <>r  the  ring  interferometer  sensitivity  A heterodyne 
technique  evolving  from  the  haste  dc  conligurculion  shown 
m log  I 


To  illustrate  the  modulation  with  homodyne  scheme  as  shown  in  Fig.  II,  we  propose  to 
use  magneto-optical  modulators  to  provide  (a)  the  dc  90°  bias  as  discussed  earlier,  and  (b)  an 
additional  modulation  around  \<t>  — tt/2  by  imposing  a sinusoidal  signal  \H  with  a proper 
amplitude  such  that  the  modulation  is  linear  in  <b.  An  alternative  to  (a)  and  (b)  is  to  use  a sin- 
gle large  ac  modulation,  producing  a phase  shift  up  to  rr , and  then  to  count  and  record  the 
maxima  in  the  output  signal.  Each  maximum  corresponds  to  the  phase  shift  being  scanned 
across  tt/2,  the  maximum  sensitivity  point.  Pursuant  to  the  desirable  use  of  near-infrared 
wavelengths,  we  point  out  that  many  promising  magneto-optical  modulators  based  on  E traday 
effect  have  been  devised  and  are  ready  for  adoption  here  155],  We  emphasize  that  only 
through  modulation,  i.e.,  ac  operation,  can  the  formula  given  by  Moss,  et  al  (6)  relating  the 
sensitivity  to  the  photon  noise  be  used.  In  the  dc  bias  case,  the  challenge  is  to  maintain  the 
dc  bias  to  the  highest  possible  stability. 
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l'ig  1 1 — Another  example  of  ac  dclcvllon  scheme:  A 
phase  lie  bias  with  modulation  or  a large  ac  modulation 
alone. 


In  addition  to  the  above-mentioned  techniques,  there  is  yet  another  one  which  is  based  on 
pulse  operation  (Fig.  12).  The  key  element  is  the  reciprocal,  gated  phase  shifter  which  is  an 
electro-optical  gate.  The  advantage  of  this  scheme  lies  in  the  fact  that  existing  integrated  optics 
techniques  can  be  used,  although  these  techniques  can  also  be  developed  for  the  first  two 
schemes. 


Kig.  12  — Another  possible  example  of  ac  detection 
scheme:  A reciprocal,  galcd-phasc  shifter,  amplitude- 
detection  technique. 
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SENSITIVITY  LIMITS 


Close  investigation  of  Fig.  6 suggests  the  following  question:  Given  a total  loss  rate  and 
optical  wavelengths,  does  there  exist  an  optimum  optical  fiber  length  to  produce  a minimum 
detectable  rotational  rate?  We  also  see  from  Eq.  (33)  that  S/N  is  proportional  to  /'cw  and,  in 
turn,  to  P(h  the  laser  output  to  the  ring  interferometer.  Increasing  P0  will  improve  the  S/N  for 
a given  L.  However,  there  exists  a limit  beyond  which  the  nonlinear  optics  effects,  such  as 
stimulated  Brillouin  scattering,  will  become  dominant.  If  no  upper  limit  is  observed,  fiber  dam- 
age will  result.  Therefore,  the  practical  sensitivity  of  a fiber  ring  interferometer  is  limited  by 
spontaneous  noises  at  low  power  and  by  stimulated  noises  at  higher  power.  Because  knowledge 
of  these  limits  is  extremely  relevant  in  the  preliminary  design  stage  of  a ring  interferometer, 
inasmuch  as  the  fiber  length  and  power  level  are  the  important  factors  that  determine  the  final 
package  size  and  unit  cost,  we  will  assess  these  limits  here. 

Optimum  Length  and  Minimum  Detectable  Rotational  Rate  at  Low  Power 


We  consider  an  improved  version  of  the  optical  arrangement,  in  that  a stable  90°  dc  phase 
bias  has  been  introduced.  With  such  an  introduction,  Eq.  (8)  for  a single  detector  becomes 


'r 


D , 


+ Ph)  ~ T Ps w cos  A*  + 


y(/»cw  + Pn>  + y P cw  sin  (A <t>) 


(34) 


where 


A0  - 4>cw  - <t>c, 


4 nRL 

Ac 


ft. 


If  the  dc  term  is  filtered  out,  the  power  of  the  signal  is 

•S’  - -7  DI 2G2P2V  sin2  (A </>)RL 

4 

and  the  quantum-noise-limited  S/N  will  be 
s yD2<72/>2,  sin2  (A^)RL 
~N  ~ *<72(Z>(/>CW  + PN) ]/?L*o 


(35) 


I _IL_ 

4 hvB0 


+ /V 


sin2  (A <t>) 


= F(L)  sin2 
or,  we  can  write 


4 nRL 


Ac 


ft 


(36) 


sin  (0ftL)  - [(S/N)  (l/F(L))V/2, 
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whence 

11  - sin  1 Hl/F(L))(S/N)V'*. 

The  optimum  fiber  length  L0  is  obtained  by  setting  (dfl/dL)Lm,  = 0 with  S/N  - I. 
Let  U = [F(L)\  l/2;  Eq.  (37)  then  becomes  (for  S/N  - I) 


a -jj-  sin-'  luan 


and 


dil  1 


dL  pL 2 

By  setting  (di\/dL)  t a — 0,  we  obtain 
and 

,-tj 


- sin  U + 


dU 

dL 


Vl  - U 2 


sin  'Ua 


U o-  l/(Lo) 

4 dU  i 

’ Vl  - t/„2  dL 


Use  of  the  distribution  laws  in  Eq.  (32)  gives 


tV 


4/t^ol'/2  K2 


T)H 


e2  [1  + (/4  + fl  + C)Ln 


+ («'  + C')(e"r'"  - e "rL"))ul 


and 


where 


dU 

dL 


- Un~-  |1  + H(L{)\, 


//(£.„)  - 


A+B+C 


+ (fl'+C)  (e"f'"  + e 


1 + (A+B+C)L„  + (fl'+C)  (e"2'”  - f "2'")  ' 


and 


^ j qsri  0 ■ 0.2«sr,  C ■ 0.2fi^ 


(37) 


(38) 


(39) 


(40) 

(41) 


B - 0.1  — . C'  - 0.1  — . 

«T  “T 

if  the  bare  cladding  fiber  is  used  as  an  example.  Appropriate  corresponding  values  must  be 
used  for  a mode-stripped  fiber,  but  the  procedure  will  remain  the  same.  Substituting  Eq.  (40) 
and  (41)  into  (39)  gives  a condition  for  the  determination  of  L0for  a given  value  of  the  param- 
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eter  fi.  Equation  (39)  can  only  be  solved  numerically.  We  have  carried  out  numerical  solu- 
tions for  the  fused  silica  fiber  discussed  above  for  both  A = 0.633  and  A - 1.1  /xm 
Numerical  values  of  the  other  variables  needed  in  the  evaluation  are  given  in  Table  2.  For  a 
fiber  coil  of  radius  R — 10  cm,  numerical  solutions  are  that  (a)  for  A — 0.633  /xm,  the 
optimum  length  Lu  — 1.56  km,  and  (b)  for  A = I I /xm,  L„  - 14.3  km.  The  corresponding 
minimum  rotational  rates  flmin  are  3.79  x 10  ® rad/s  for  A - 0.633  /xm  and  1.2  x 10  8 rad/s 
for  A “ 1.1  /xm.  These  are  listed  in  Table  3. 

As  we  expected,  with  due  consideration  of  light  scattering  in  fiber  and  90°  phase  bias 
improvement,  our  results  indicate  less  reason  for  optimism  than  found  by  in  Vali  and  Shorthill 
[3).  However,  we  think  that,  with  the  implementation  of  any  of  the  ac  detection  schemes  and 
heavy  mode  stripping  mentioned  above,  we  may  achieve  an  order  of  magnitude  improvement 
in  sensitivity. 


Table  2 — Numerical  Values  of  the  Parameters  used 
to  Evaluate  the  Optimum  Fiber  Length  (An  Example) 


Case 

fl 

flo 

(Hz) 

Pc 

(mW) 

“i 

C 

n 

(W) 

m 

nmw 

0.5 

1 

2 

0.5 

0.5 

irrrai 

ign'irm 

0.1 

1 

2 

0.5 

0.5 

firmi'M 

HHTM 

Table  3 — Quantum-Noise-Limited  Performance* 


A 

(/xm) 


“r 

(dB/km) 


L0 

(km) 

Low-Power 


Minimum  Detectable  0 
(deg/h) 


Low-Power 


High-Power 


0.633 


1.56 


(P„-  2 mW) 
0.0078 


</>„  - 81.2  mW) 
0.0009 


1.1 


0.43 


14.3 


(P„  - 2 mW) 
0.0025 


(P„  - 14.4  mW) 
0.0007 


*90”  phase  bias  and  liber  scattering;  K - 10  cm 


Effect  of  Stimulated  Brillouin  Scattering  at  High  Power 

Previous  works  (1,3,4)  dealing  with  optical  fiber  ring  interferometers  and  our  discussions 
so  far  are  limited  to  linear  optics.  However,  the  S/N  will  improve  as  we  increase  the  laser 
power,  hence  the  power  coupled  into  the  fiber,  as  evidenced  by  Eq.  (36).  The  sensitivity  in 
terms  of  the  detectable  rotational  rate  will  also  improve  because  1 

sin  '£/„«  Un  for  small  (/<>  But  U0  « II  1/2  * P„  l/2,  hence  ftmm  « P0  1/2  The  upper  limit  for 
P„  will  be  at  a level  where  noises  induced  by  stimulated  processes  become  dominating  and  start 
to  degrade  the  S/N  and  ft  min.  These  nonlinear  optics  effects  will  become  more  important  as  we 
reduce  the  total  losses  of  optical  fiber,  for  then  more  power  input  to  these  fibers  would  be  avail- 
able for  stimulation. 
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The  stimulated  Raman  and  Brillouin  scattering  in  an  optical  fiber  have  been  investigated 
theoretically  by  R.G.  Smith  (561.  Based  on  the  assumption  of  nondepleting  pump  beam  and 
the  criterion  that  the  critical  power  is  reached  when  the  scattered  power  in  the  stimulated  pro- 
cess is  equal  to  the  local  pump  power.  Smith  establishes  that  the  backscattcring  stimulated  Bril- 
louin process  requires  the  least  pump  power.  Therefore,  we  shall  concentrate  on  this  process. 


In  order  to  present  a discussion  appropriate  to  the  operation  of  fiber  ring  interferometers, 
we  adopt  a new  criterion  for  the  critical  power.  It  stales  that,  for  a chosen  fiber  length,  the  critical 
power  is  the  input  laser  pump  power  for  which  the  backscattered  Stokes  wave  acquires  a power  level 
as  high  as  the  transmitted  power  of  the  beam  coming  from  the  other  end  of  the  fiber.  Let  Pn  (0)  be 
the  pump  power  at  the  end  x-0,  P„(L)  be  that  at  x-L,  P,( 0),  the  backscattered  Stokes  wave 
emerging  from  x— 0.  Then,  mathematically,  the  criterion  translates  into 

#»,(' 0)  - Pn(L)e'nTL.  (42) 

Since  we  have  assumed  symmetry  in  the  CW  and  CCW  directions,  we  can  set  PU  \ - P,,( 0) 
such  that 


P,(0)  - P„(0)  e " T>\  (43) 

According  to  Smith  (56),  the  backscattered  stimulated  Brillouin  scattering  has  the  lowest  thres- 
hold and  can  be  shown  to  follow  the  formula 


P,(0) 


■Jn 

2 


( kT ) (Av«)  -Jy-P-^~a  7.' 


x 


—2a  + 


>0^,(0) 

Aar 


(44) 


where  v,  and  v„  are  the  frequencies  of  the  Stokes  wave  and  the  acoustic  phonon,  respectively, 
H is  the  line  width  of  the  spontaneous  Brillouin  scattering,  T is  the  absolute  temperature  of 
the  fiber,  y0  is  the  peak  gain  coefficient  for  the  stimulated  Brillouin  process,  A is  the  core 
cross-sectional  area.  This  expression  is  applicable  for  all  a , Ly  since  we  have  used  Smith’s 
expression  for  full  gain  given  by  his  Eq.  (B5).  Substitution  of  (44)  into  Eq.  (43)  gives  a final 
equation  for  evaluating  the  critical  input  pump  power  when  all  the  other  parameters  are  known. 


Two  cases  have  been  evaluated:  namely,  X - 0.633/um  and  X - 1.1  /*m.  In  doing  so,  we 
have  ignored  the  frequency  shift  of  the  scattered  light  and  set  vs  - i/la,cr.  We  also  arbitrarily 
take  L — Lq,  which  are  determined  from  low-power  calculations.  Other  values  used  in  the 
evaluation  are  listed  in  Table  4. 


Results  show  that,  for  the  case  X - 0.633  /un,  the  critical  input  pump  power  Pr( 0)  - 10 
mW,  and  for  X «•  1.1  /im,  PP( 0)  — 1.8  mW.  In  terms  of  the  laser  output  power 
P p(0)/ll/2(l— at)C]),  we  have  P0  - 1.8  mW  x 8 - 14.4  mW  for  X -1.1  gm  and 
/*0  — 81.2  mW  for  X — 0.633  *im. 

Based  on  our  criterion,  the  quantum-limited  S/N  ratio,  according  to  Eq.  (36),  with 
PN  - Pcw,  is  given  by 

1^1  1 V^CW  . 
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where  we  have  set  PJ 0) 
rate  is  given  by 


P\  ” PLW-  Setting  S/N  = 1,  the  minimum  delectable  rotational 


(l 


min 


vP„(oi 


x 

a 


(46) 


where  the  approximation  sin  'x~x  for  small  x,  and  /'JO)  = /\w  = P,,(0)c  is  assumed.  As 
an  example,  let  us  consider  the  same  ring  interferometers  used  for  low  power  rates  as  above. 
Our  numerical  calculations  show  that,  for  the  case  X = 0.633  /or,  llmin  — 4.44  x 10  9 rad/s, 
or  9.16  x 10  4 deg/h,  whereas  for  X =-  1.1  flm,n  - 3.35  x 10  ’ rad/s,  or,  6.9  x 10  4 
deg/h.  These  values  are  listed  in  the  last  column  of  Table  3.  it  is  of  interest  to  compare  these 
with  the  counterparts  for  the  low-power  cases:  For  X — 0.633  fim,  the  high-power  performance 
in  rotational  rate  sensing  has  about  88%  improvement  over  the  low-power  value.  For 
X — 1.1  pm,  we  gain  about  72%  imnrovement.  The  optimum  fiber  length  L„,  which  gives  flm,„ 
for  S/N  * 1,  can  be  determined  from  Eq.  (46).  The  result  is  that  L0  - 2/a  ,.  This  condition 
still  applies  if  we  neglect  all  the  scattered  light  components  (PN  = 0),  and  it  is  different  from 
the  condition  L()  - 0.87/a,  as  given  by  Ref.  3. 


What  is  the  sensitivity  for  the  ideal  case?  This  question  may  be  become  relevant  when  the 
state  of  fiber  development  reaches  a high  level  of  perfection  whereby  the  scattered  light  is 
insignificant  and  when  ultrastable  dc  bias  and  ac  modulation  techniques  can  be  implemented. 


In  the  absence  of  scattered  light,  the  ideal  S/N  for  the  dc  detection  scheme  is 


whence  (for  S/N  - 1) 


5 

N 


i 
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vPc» 


hvB{) 


sin2  (/3 LSI) 


(47) 
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min 


4 hvB0 
VP„(0) 


1/2 


(48) 


Comparing  with  (46),  we  see  that  the  minimum  detectable  rotational  rate  in  the  ideal  case  is  a 
factor  I/V2  lower  than  the  high-power  case.  If  the  differential  scheme  is  used,  the  ideal 
minimum  detectable  rotational  rate  is  a factor  of  one-half  lower.  For  the  heterodyne  case,  we 
would  expect  a factor  of  l/(2\/2  )lower. 


CONCLUDING  REMARKS 

Optical-fiber,  ring  interferometers  and  their  variants  can  also  be  used  for  the  detection  of 
other  physical  parameters  besides  being  highly  sensitive  rotational  rate  transducers.  Physical 
effects  such  as  the  acoustic  pressure  field,  stress  distributions  in  composite  materials,  the  tem- 
perature field,  and  magnetic  fields  can  be  detected  through  the  optomechanical,  magneto-optical 
properties  of  the  fiber  materials.  Acoustic  pressure  field  measurements  have  been  demon- 
strated [51,  and  other  effects  are  currently  under  evaluation  in  our  laboratory.  The  present  sen- 
sitivity analysis  will  certainly  facilitate  the  evaluation. 

Up  to  now,  our  discussions  have  not  included  the  depolarization  effects  in  a single-mode 
fiber.  In  certain  applications,  such  as  stress  measurements  or  temperature  measurements,  the 
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polarization  state  of  the  optical  beam  in  the  fiber  can  be  used  to  our  advantage;  but,  in  ring 
interferometers  for  rotational  rate  measurements,  if  a linearly  polarized  beam  is  used,  depolari- 
zation will  cause  signal  level  reduction,  (a  factor  of  2 for  complete  depolarization).  Therefore, 
the  effect  of  depolarization  due  to  uncontrolled  nonrecriprocal  changes  in  the  fiber  should  be 
investigated.  Nonlinear  effects  with  pump  depletion  should  also  be  investigated.* 

Our  analysis  of  the  sensitivity  has  shown  that  when  laser  input  power  is  low,  spontaneous 
noises  are  important,  and  they  degrade  the  S/N  ratio.  The  higher  the  fiber  loss,  the  stronger 
the  scattered  light,  hence  the  noise  power.  However,  the  level  of  degradation  may  be  reduced 
considerably  if  a proper  mode  stripper  can  successfully  remove  these  cladding-guided, scattered 
components.  We  can  expect  two  orders  of  magnitude  in  improvement.  But,  when  the  input 
power  is  increased  to  above  certain  levels,  the  nonlinear,  stimulated  scattering  will  dominate. 
The  power  level  beyond  which  nonlinear  optics  effects  are  important  is  lower  when  optical  fiber 
of  lower  loss  is  used.  On  the  other  hand,  sensitivity  will  increase  as  the  square  root  of  power. 
Our  analysis  shows  that,  when  the  stimulated  Brillouin  scattering  is  taken  as  the  sole  source  of 
noise,  the  gyro  sensitivity  is  better  than  that  of  the  low-power  case.  We  conclude  that,  in  pur- 
suant to  the  low-loss,  single-mode  fiber  development,  the  fiber-ring  interferometer  should 
operate  with  the  highest  possible  input  power. 

Because  the  optical  fiber  technology,  light  sources  such  as  laser  diodes  and  detectors  are 
currently  under  intense  research  and  development  for  optical  fiber  communications  systems,* 
the  promising  use  of  optical  fiber  as  the  transducer  can,  in  fact,  be  tied  to  the  advancement  in 
that  field.  Researchers  in  optical  fiber  communications  systems  have  pointed  out  the  clear 
advantage  of  operating  the  systems  in  the  infrared  [49,56).  From  the  point  of  view  of  scatter- 
ing noise,  ring  interferometers  should  also  follow  this  trend.  What  seems  to  be  the  limitation 
today  may  soon  be  removed  tomorrow.  Our  optimistic  viewpoint  is  supported  by  the  favorable 
comparison  of  the  ring  interferometer  with  a specific,  well-developed  mechanical  gyro  and  an 
advanced  ring  laser  gyro.  This  is  shown  in  Table  5.f  We  conclude  that  ring  interferometers  as 
fiber  gyros  can  soon  be  the  alternative  to  ring-laser  gyros  in  certain  applications.  The  ultimate 
sensitivity  as  a gyro  is  anticipated  to  be  even  higher  than  the  ring  laser  gyro. 
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